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Abstract The origins of the hot solar corona and the supersonically expanding solar wind 
are still the subject of much debate. This paper summarizes some of the essential ingredients 
of realistic and self-consistent models of solar wind acceleration. It also outlines the major 
v^ ' issues in the recent debate over what physical processes dominate the mass, momentum, 

and energy balance in the accelerating wind. A key obstacle in the way of producing real- 
istic simulations of the Sun-heliosphere system is the lack of a physically motivated way 
of specifying the coronal heating rate. Recent models that assume the energy comes from 
C^ ■ Alfven waves that are partially reflected, and then dissipated by magnetohydrodynamic tur- 

(-H I bulence, have been found to reproduce many of the observed features of the solar wind. This 

Oj, paper discusses results from these models, including detailed comparisons with measured 

' ' plasma properties as a function of solar wind speed. Some suggestions are also given for fu- 

t~ , ture work that could answer the many remaining questions about coronal heating and solar 

wind acceleration. 



Keywords solar corona • solar wind • turbulence • waves 



1 Introduction 



> 

Qs , It has been known since the early part of the twentieth century that the temperature in the 

(^ • Sun's outer atmosphere undergoes a rapid inversion, from a relatively cool {T < 10^ K) pho- 

tosphere and chromosphere to a hot and ionized (T > 10^ K) corona. Since the dawn of the 
f-^ ■ space age, coronal heating has also been known to be linked to the acceleration of a plasma 

f^ I outflow that reaches speeds of 300-800 km s^' in interplanetary space. In the years since 

the twin problems of coronal heating and solar wind acceleration were formulated, many 
physical processes have been suggested to be responsible. Only a small fraction of the me- 
chanical energy in the Sun's sub-photospheric convection zone needs to be converted to heat 
p\^ • in order to power the corona. However, it has proved exceedingly difficult to distinguish be- 

J-j , tween competing theoretical models using existing observations. Recent summaries of these 
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iroblems and cont r oversi e s have been pre sente d by, e.g ., AschwanderJ ll2006h . iKlimchukI 



2006h . lZurbuc"to] bOOTh . lHollw^ ( l2008h . and lCranmej (|200^ . 

Most of the proposed physical models can be grouped into two broad paradigms. As will 
be seen below, the basic difference between these paradigms concerns the overall topological 
connectivity of the magnetic flux tubes that feed the solar wind: 

1 . If solar wind flux tubes are open to interplanetary space — and if they remain open on 
timescales comparable to the time it takes plasma to accelerate into the corona — then 
the main sources of energy must be injected at the footpoints of the flux tubes. Thus, 
in wave/turbulence-driven (WTD) models, the convection-driven jostling of the flux- 
tube footpoints is assumed to generate wave-like fluctuations that propagate up into 
the extended corona. These waves (usually Alfven waves) are often proposed to par- 
tially reflect back down toward the Sun, develop into strong magnetohydrodynamic 
(MHD) turbulence, and dissipate gradually. These models have been shown to natu- 
rally produce realistic fast and slow wind conditions with wave amplitu des of the same 
order of m agnitu de as those observed in the corona and heliosphere I HoUwed 19861: 



VellietalJll9 91: Wang and Sheelev 1991; Matthaeus et al. 1999: Suzuk i and Inutsuk3 



2006 : ICranmer et al..2007. : ,Wang et al..2009, : ,Verdini et al..20ia : ,Matsumoto and Shibatal 



201(3) 



2. Near the Sun, all open magnetic flux tubes are observed to exist in the vicinity of 
closed loops that a re evolving on a wide ran ge of spatial and time scales in a complex 
"magnetic carpet" ( ITitle and Schriiveiill998h . Thus, it is natural to propose a class of 
reconnection/loop-opening (RLO) models, in which solar wind flux tubes are assumed 
to be influenced by magnetic reconnection with closed-field regions that are continu- 
ously emerging, fragmenting, and being otherwise jostled by convection. In these mod- 
els the mass, momentum, and energy of the solar wind is input from loops of varying 
properties in the low corona. Some have sugge sted that RLO-type energy ex c hange pri- 
marily occurs on small, supergranular scales JAxford and McKenzielll992l : iFisk et al.l 
ll999l : lFiskll2003l : ISchwadron and McComasll2003h . However, other models have been 
pro posed in which the reconnection occurs in and between large-scale coronal stream- 
ers l lEinaudi et aO 19991 : ISuess and Nemevll2004l : lAntiochos et al.ll2010h . 

In the interest of brevity, this paper does not review the wider range of empirically based 
solar wind models that do not contain self-consistent coronal heating p hysics. These mod els 
are often quite sophisticated in their treatment of multi-fluid (e.g., Ha nsteen et al. 19971: 



2003; Lie-Svendsen and Esser 2005) or multi-dimensional (e.g., Endev e et al.ll2003L l2004l : 



Roussev et al. 2003 ; Rilev et al. 2006 ; Nakamizo et al. 2009) effects. However, their use of 



ad hoc rates of heating and momentum deposition places them in a separate category from 
the self-consistent models discussed here. 



2 Essential Ingredients 

Before reviewing the results of WTD or RLO models, it is useful to summarize what kinds 
of "physics inputs" are necessary to build a self-consistent model of solar wind acceleration. 
The first models of the solar wind were not self-cons istent, but the y contained many 
valuable insights that led the way to future improvements. IParkerl (Il958h found that a con- 
stant coronal temperature of order 10*" K (i.e., roughly consistent with strong radial electron 
conduction) provides enough of an outward gas pressure gradient force to overcome grav- 
ity and produce a natural transition from a subsonic (bound, negative energy) state near 



the Sun to a supersonic (escaping, positive energy) state in interplanetary space. Through- 
out the 1960s, new models with different prescribed coronal temperatures T{r) explored 
the parameter space of possible time-steady solar wind solutions. Some of these mod- 
els included tempe r atures consistent with polytropic equations of state (i.e., P ex pi'), and 
iHolzer and AxfordI 1 1970l ) found that y < 1.5 is necessary for an accelerating solar wind. 



Also. lSturrock and Hartld ( 1 196a) solved two-fluid (Tp ^ T^) energy equations with heat con- 
duction, and found that some kind of "extra" energy addition is needed in the corona to heat 
the plasma to the temperatures seen at 1 AU. 

The order-of-magnitude amount of coronal heating that is needed to produce the solar 
wind can be estimated from an approximate version of the internal energy conservation 
equation (see, e.g.. iLeer et al.iil982) . The time-steady version of this equation. 



, ..' 5P GMc. 

Fheat + Fcond + PU -T- + 7i 

' 2 2p r 



Grad , (1) 



describes the balance between an imposed heat flux (Fheat)> conduction along the magnetic 
field (Fcond)> radiative losses (2rad)> and fluxes of kinetic energy, enthalpy, and gravitational 
potential energy (the three terms in parentheses). When studying how the energy budget 
varies from the corona to 1 A U, the dominant terms ar e the imposed heating, the kinetic 
energy flux, and gravity. Thus. lHansteen and Leen ( Il995r) found that one can estimate 

^heat = iFheatI « (P«)corona (^ + y) (2) 

(see also iLeer and Marschlll999l : ISchwadron and McComasll2003h . Thus, if we can spec- 
ify the mass flux in the corona, pu, the escape velocity from the solar surface, Vesc = 
{IGMq/RqY'^, and the asymptotic, or terminal outflow speed, u„, we can estimate how 
much heat must be deposited in the corona. Using typical values d Cranmer et alj|2007h for 
the fast solar wind associated with coronal holes, Fheat » 8 x 10^ erg cm^^ s^' . For the slow 
solar wind associated with streamers and active regions, Fheat « 3 x 10^ erg cm^^ s^'. 

In order to go beyond order-of-magnitude estimates, models must begin to include self- 
consistent descriptions of specific physical processes involved in the mass, momentum, and 
energy conservation of accelerating plasma parcels. There are at least five essential ingredi- 
ents to such a self-consistent picture of the solar wind: 

1 . Physically motivated coronal heating. The actual origin of the imposed heat flux Fheat 
(or the equivalent volumetric rate Qheat = |V -Fheatl) must be included explicitly. The key 
difference between the WTD and RLO paradigms, as described in the previous section, 
rests on whether the heating is deposited by fluctuations within an open flux tube or via 
impulsive injection from other surrounding flux tubes. Models can differ, of course, in 
the level of detail given in the self-consistent heating rates. Some descriptions of turbu- 
lence and magnetic reconnection specify only the energy that is input into the system on 
the largest spatial scales. These models thus employ various phenomenological assump- 
tions about how that energy is eventually dissipated. Other models follow the detailed 
microphysics of the dissipation itself — usually by simulating either particle-particle col- 
lisions or wave-particle interactions. 

2. Additional momentum sources in the fast wind. In the 1970s and 1980s, it became in- 
creasingly evident that the maximum mean plasma temperature in the open-field corona 
(i.e., approximately {Tp + Ff )/2) does not exceed ~ 2 x 10*' K. Even the measurements 
of higher values of Tp in coronal holes made in the 1990s did not substantially change 



this constraint on the one-fluid average (see, e.g.. lKohl et alJcOOd : ICranmenl2009l) . A 
solar wind with this mean temperature that is driven only by ga s pressure cannot accel - 
erate to the highest speeds (700-800 km s^') measured at 1 AU ( ILeer and Holzenl980l) . 
MHD fluctuations have been shown to exert an additional "wave pressure" or outward 
ponde romotive force on the mean fluid ( lDewanll970l : lBelchenll97ll : IOfman and Davilal 
Il998r) . This appears to be a necessary component of fast-wind models, and it may even 
provide as much as half the acceleration in flux tubes connected to polar coronal holes 
(e.g., ICranmer 2004) . 

3. A self-regulating mass flux. The Sun's mass loss rate, M w 2 x 10^'^ M0 yr^^ is gen- 
erally believed to be determined in the transition region by a balance between do wn- 
ward heat conduction, local radiative losses, and the upward e nthalpy flux (Ham men 
1 19821 : IWithbro3l 19881 : iHansteen and Lee3l 19951 : ITeer et alj|l998h . Self-consistent mod- 
els should not artificially fix the properties of the transition region and low corona. In- 
stead, models should allow these properties to "float" until a natural, stable, and time- 
steady solution for the energy balance is found. 

4. Extended conduction and heating to 1 AU. Even far above the corona, in situ measure- 
ment s show that ene rgy deposition and conductive energy transfer are still occurring 
(e.g.,'Coleman 1968). The radial gradients of the proton and electron temperatures are 
substa ntially shallower than would be the case if plasma parcels were expanding adiabat- 
ically JMarsch et al.l 19831 . 1 1989l : lRichardson et al.ll995h . Helios measurements of radial 
growt h of the proton magnetic moment between 0.3 and 1 AU JSchwartz and MarscH 
Il983r) point to specific coUisionless processes that continue to affe ct the energy budget 
far from the Sun. R ecent empirical studies of the fast solar wind JBreech et al.ll2009l : 
ICranmer et al.ll2009l) have shown that it is important to take careful account of electron 
heat conduction in order to determine how the heating is partitioned between protons 
and electrons. Self-consistent models should be able to describe both the extended heat- 
ing that is observed and the transition from coUisional to coUisionless conduction in the 
heliosphere. 

5. Funnel-type magnetic field expansion. Most of the plasma that eventually is accelerated 
outwards as the solar wind seems to originate in the lanes and vertices between su- 
pergranular network cells in the chromosphere. As height increases, the strong vertical 
magnetic field decreases and the flux tubes expand laterally. Thus, the flux tubes become 
magnetic funnels that even tually merge with one another into a topologically complex 
"cano py" in the low corona ( lGabrielll976l : lDowdv et alJl98q : ICranmer and van BallegooiienI 
|2005|) . Several semi-empirical studies of coronal heating and solar wind acceleration 
show that this kind of funnel-like flux tube expans ion is necessary to producing realistic 
emission-line spectra in the transition region (e.g.. lEsser et al.l2005l : lMarsch et al.l2008l : 
iBvhring et al.ll2008l : |Pucci et al.ll2010l) . 



The ways in which the above ingredients interact with one another to produce a time-steady 
solar wind are complex and nonlinear. Even though the wind flows upwards, sometimes the 
physics of the extended corona (e.g., heating at heights of 1-2 solar radii above the surface) 
can have a significant feedback on lower regions in the atmosphere (e.g., frozen-in charge 
states that are set just above the transition region). It is important for a model to allow 
these various pieces of physics to evolve together toward a stable steady state and not be 
constrained by input assumptions. For example, the lower boundary of the model should 
not be so high as to exclude the transition region and upper chromosphere, and the upper 
boundary should not be so low as to exclude the sonic or Alfvenic critical points of the flow. 



3 Wave/Turbulence Models 

There has been substantial work over the last few decades devoted to exploring the idea 
that coronal heating and solar wind acceleration may be explained as a result of the dissipa- 
tion of waves and turbulent fluctuations. No matter the relative importance of reconnections 
and loop-openings in the low corona, we do know that waves and turbulent motions are 
pr esent eve rywhere fro m the photosphere to th e heliosphere (see observational summaries 
of lXu and Marsch 1995uCrarmier and van Ballegooiien 2005 : Aschwanden 2008). Thus, it 
is of interest to determine how waves affect the mean state of the plasma in the absence of 
any other sources of energy. 

ICranmer et al. (2007) described a set of WTD-type models in which the time-steady 
plasma properties along a one-dimensional solar wind flux tube are computed. These model 
flux tubes are rooted in the optically thick solar photosphere and are extended into interplan- 
etary space. The numerical code developed in that work, called ZEPHYR, solves the one- 
fluid equations of mass, momentum, and energy conservation simultaneously with transport 
equations for Alfven and acoustic wave energy fluxes. ZEPHYR is the first code capable of 
producing self-consistent solutions for the photosphere, chromosphere, corona, and helio- 
sphere that combine: (1) shock heating driven by an empirically guided spectrum of acous- 
tic waves, (2) extended heating from Alfven waves that get partially reflected and then are 
dissipated by MHD turbulence, and (3) wind acceleration from gradients of gas pressure, 
acoustic wave pressure, and Alfven wave pres sure. 

The only input "free parameters" to the Cranmer et al.l ( 120070 models were the photo- 
spheric lower boundary conditions for the waves and the radial dependence of the magnetic 
field along the flux tube. Photospheric measurements of the horizontal motions of intergran- 
ular flux concentrations (i.e., G-band bright points having S ^ 1 .5 kG) were u sed to constrain 
the fr equency spectrum of Alfven waves at the lower boundary (see also iNisenson et al.l 
l2003h . All models shown below used the same lower boundary conditions and differed only 
in the rates of flux-tube superradial expansion. 

The self-consistent coronal heating in the ZEPHYR models is the result of propagat- 
ing Alfven waves being partially reflected by radial gradients in the density and magnetic 
field strength. It has been shown that once there are counter-propagating wave packets 
that interact with one another along a flux tube, a nonlinear turbulent cascade can then 
occur relatively quickly (Irosh nikoMll964l : iKraichnan 1965) . The energy flux in the cas- 
cade, from large to small eddies, terminates in dissipation and heating that extends from 
the low corona out into the heliosphere (see also Matthaeus et al. 1999, Dmitruk et al. 2002; 
Cranmer and van Ballegooiien 2005; Chandran and HoUweg 2009; Verdini et al. 2009, 2010; 
CranmeJbOld) . 



In addition to the properties of the waves and turbulence, there are other aspects of the 
models that determine how much coronal heating occurs. One of these is the radial location 
of the Parker "critical point." This is the point at which the wind speed exceed s a criti- 
cal sp eed defined by the sound speed and the MHD wave amplitudes (see, e.g., IjacquesI 
Il977h . In models where the magnetic flux tubes expand purely radially, there is usually 
just one unique location for this critical point. However, in flux tubes that undergo super- 
radial expansion there are multiple possible locations where the critical point could be lo- 
cated. These lo cations correspond to loc al m inima in a potential-e nergy-like quantity that 
was defined by Koop and Holzen (Il97a) and I Vasquez et al.l (l2003h . However, only one of 
these points corresponds to a global energy minimum, and thus only this one critical point 
location gives a stable and time-steady solar wind. For some models, changing the flux- 
tube expansion only slightly is enough to alter the relative depths of these potential en- 



ergy wells, such that the critical point shifts abruptly from a location close to the Sun 
to one much farther from the Sun. When this occurs the global mome ntum and energy 
balance of the solar wind changes abruptly as well. Early studies (e.g., iLeer and Holzen 
1 1 980l : IPneumarJI 1 980n showed that high critical points — where most of the heating is in the 
subsonic low corona — correspond to dense and slow solar winds. Conversely, low critical 
points — where the heating is mainly in the supersonic outer corona — correspond to low- 
density a nd fast solar winds. Fo r a "stretched dipole" model of the solar-minimum magnetic 
field, the ICranmer et al.l ( 120071) models show precisely this kind of discontinuity between 
fast and slow winds a t a heliospheric latitude of ~ 20°, similar to what Ulysses observed 
(lGoldst einetal.!l l996h . 

Figure[T]s hows several different comparisons between solar wind measurements and the 
results of the ICranmer et al.l ( 120071) models. In all six panels the observations and models 
are sorted by the solar wind speed at 1 AU. In panel (a), in situ proton temperature data 
are taken from t he ACi^/SWEPA M online archive for the time periods between 1998 and 
2005 (see also Matthaeus et al.ll2006h . Note that the plot juxtaposes the measured Tp val- 
ues with the ZEPHYR one-fluid mean temperatures {{Tp + T^)/T), so th e comparison is not 
exact. In panel (b), the coronal superradial expansion factor defined by IWang and SheelevI 
1 1 1 9901) is shown as a func tion of wind speed, along with the empirical trend of anticorrela- 
tion also demonstrated bv IWang and SheelevI ( 11990') . In panel (c), ACE proton temperatures 
and densities have been combined to form the specific entropy quantity that was found by 
IPagel et al.l 1 12004 1 to correlate strongly with wind speed. In pa nel (d), th e power in magnetic 



fluctuations measured by Helios between 0.3 and 0.5 AU bv iTu et al.l ( |I992 | ") is compared 



, tne po 



with a similar quantity estimated from the models in the way described by ICranmer et al.l 
(l2007h . 

Panels (e) and (f) of Figure[T]show ion data from UlysseslSi^lCS taken at two different 
phases of that mission (light gray: 1990-1994, dark gray: 1994-1995). In panel (e), the 
measured ratio of 0^+ to 0*+ number density is compared to ZEPHYR models of this 
"frozen in" nonequilibrium ionization state. Although the overall trend with wind speed is 
similar to that in the data, there is an overall shift downward in the modeled ratios. This may 
be due to the fact that the models assume the electron velocity distribution to be Maxwellian; 
i.e., the models ignore the a dditional ionization th at would be caused by a suprathermal 
"electron halo" in the corona dEsser and Edgan l2000). In panel (f), the measured ratio of iron 
to oxygen elemental a bundances — no rmalized to their photospheric values — is compared to 
models that apply the lLamind ( l2004h idea for first-ionization-potential (FIP) fractionation. 
This idea utilizes the ponderomotive force exerted by Alfven waves to accelerate ions in 
the partially ionized upper chromosphere, while leaving the neutrals unaffected. Note that 
the model results shown in panels (e) and (f) contradict the commonly held assertion that 
slow-wind FIP and charge-state properties ca n only be explaine d by the injection of plasma 
from closed-field regions on the Sun (see also lPucci et alj|2010n . 

Three other successful predictions of the ZEPHYR models are summarized below: 

1 . Recent HinodefElS measurements of strong Doppler shifts at the edges of active regions 
indicate outflow speeds of order 100 km s~^ in the coronal source regions of some 
sl ow wind s t reams (Harra et al. 2008 : iSubramanian et al.ll2010l) . As shown in Figure 7a 
of lCranmea ( I2010I) . the ZEPHYR model of an active-region-associated slow wind flux 
tube naturally exhibits a local maximum in the outflow speed of the observed order of 
magnitude a t heights between 0.02 and 0. 1 Rq above the photosphere. 

2. The original ICranmer et alj ( 120070 model of the fast wind associated with polar coro- 
nal holes used a magnetic field model consistent with the 1996-1997 solar minimum. 
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Fig. 1 Comparisons between the ZEPHYR models (thic k black cu rves) and observational data in the solar 
wind (gray regions). See the main text for details (see also lCranmer et al. 2007; Cranmer 2009). In all panels 
except (d), the gray regions show the data binned by solar wind speed, showing only the regions within ±1 
standard deviation of the mean value at each speed. 



However, the more r ecent 2008-2009 minimum has been seen to be quite different. 
ICranmer et alj ( 120 101) produced a new ZEPHYR model of the fast polar wind that used 
a lower magnetic field consistent with both solar-disk and in situ measurements taken 
during 2008-2009. The model produced changes in the plasma properties at 1 AU that 
agree well with Ulysses measurements (e.g., McComas et al. 2008). For example, in 
both the models and the measurements, the wind speed u remains relatively unchanged, 
but the density n and temperature T decrease by factors of order 20% and 10%, re- 



spectively. The decreases in gas pressure (proportional to riT) and dynamic pressure 
(proportional to nu^) are between 20% and 30% for both the observations and models. 
The Helios probes measured Faraday rotation fluctuations (FRFs) of polarized radio sig- 
nals that passed through the solar corona at heliocentric impact parameters between 2 
and 15 R,^, in the ecliptic plane. The magnitude of these fluctuations depends not only 
on the amplitude of the Alf ven waves in the cor ona, but also on the density and the 
turbulent correlation length. iHoUweg et al.1 ( 1201 0) compared the measured FRFs with 
predicted values from ZEPHYR, and found excellent agreement when using the equato- 
rial streamer model originating at a colatitude of 28° . 



4 Reconnection/Loop-Opening Models 

It is clear from observations of the Sun's highly dynamical "magnetic carpet" that much 
of the coronal heating in closed-field regions is driven by the interplay between the emer- 
gence, separation, merging, and cancellation of small-scale flux tubes. Magnetic reconnec- 
tion seem s to be the most likely ch annel for the built-up magnetic energy to be converted to 
heat (e.g. JPriest and Forbesll2000h . Thus, the RLO id ea has a natural app eal since all open 
flux tubes are rooted in the vicinity of closed loops JDowdv et al.lll98q) . In fact, isolated 
RLO-like reconn ection events are a l ready observed in co ronal holes as polar jets by SOHO 
and Hinode (e.g.. lWang et al.lll998l : IShimoio et al.ll2007l) . Also, there are observed correla- 
tions between the length s of coronal loops, the electron temperature in the low corona, and 
the wind speed at 1 AU JGloeckler et alj|2003r) that are highly sugge stive of a net transfer 
of magnetic energy from the loops to the open-field regions (see also lFisk et al.lll999l : lFislJ 
|20P3). 

Testing the RLO idea using theoretical models seems to be more difficult than testing 
the WTD idea because of the complex multi-scale nature of the relevant magnetic topology. 
It could be argued that one needs to create three-dimensional and time-dependent models of 
the magnetic carpet in order to fully take account of all interactions between the closed and 
open flux systems. Several key questions remain to be answered. For example, how much 
of the magnetic energy that is liberated by reconnection goes into simply reconfiguring the 
closed fields, and how much goes into changing closed fields into open fields? Specifically, 
what is the actual rate at which magnetic flux opens up in the magnetic carpet? Can the 
observed polar jets provide enough energy to drive a significant fraction of the solar wind? 
Lastly, how is the reconnection energy distributed into various forms (e.g., bulk kinetic en- 
ergy, thermal energy, waves, or energetic particles) that can each affect the accelerating wind 
in different ways? 

Recent work has started to provi de tentative answers to some of the above questions. 
ICranmer and van Ballegooiienl 1120 lOf) developed Monte Carlo simulations of the time-varying 
magnetic carpet and its connection to the large-scale coronal field. These models were con- 
structed for a range of different magnetic flux imbalance ratios — i.e., for both quiet regions 
and coronal holes. The models agree with observed emergence rates, surface flux densi- 
ties, and number distributions of magnetic elements. Despite having no imposed super- 
granular motions in the models, a realistic network of magnetic "funnels" appeared sponta- 
neously as the result of diffusion-limit ed aggregation from smaller magnetic concentrations. 
ICranmer and van Ballegooiienl ( 1201(1) computed the rates at which closed field lines open 
up (i.e., the recycling times for open flux), and they estimated the energy fluxes released in 
reconnection events that involve the opening up of closed flux tubes. For quiet regions and 
mixed-polarity coronal holes, these energy fluxes were found to be significantly smaller than 



those required to accelerate the solar wind. In other words, only a tiny fraction of the Poynt- 
ing flux delivered into the corona by emerging bipoles seems to be released via magnetic 
reconnection in RLO-like events. On the other hand, for the most imbalanced (i.e., unipo- 
lar) coronal holes, the energy in flux-opening events may be large enough to power the solar 
wind. However, in those cases the overall recycling times are far longer than the time it takes 
the solar wind to accelerate up into the low corona. Thus, RLO processes on supergranular 
scales may be responsible for the intermittent jets in coronal holes, but probably not for the 
majority of the "bulk" solar wind acceleration. 



5 Conclusions 

Despite recent progress made with both the WTD and RLO approaches to plasma heating 
and acceleration, we still do not have conclusive answers about whether one idea or the 
other is dominant in the actual corona and solar wind. It is possible, of course, that quali- 
tatively different mechanisms may govern the heating and acceleration in different types of 
solar wind streams. It is also possible that some kind of combination of the WTD and RLO 
paradigms may be more valid than either idea in isolation. 

An important next step in the process is to determine what specific observations can 
best test these ideas, with the goal of convincingly verifying and/or falsifying them. Also, a 
related future step will be to build three-dimensional models of the Sun-heliosphere system 
that include WTD and RLO physics as their "coronal heating functions." These kinds of 
simulations can be customized for specific time periods, and be used to make straightfor- 
ward comparisons with various kinds of existing remote-sensing and in situ measurements. 
Some recent progress in producing computationally efficient approximations to the rates of 
WTD wave refl ection and heating has been reported bv lChandran and HoUwea ll2009h and 
ICranmed ( l2O10h . 

In addition to expanding the scope of the models, it will also be important to develop a 
better understanding of the physics of MHD wave generation, propagation, and dissipation. 
Recent observations of Alfven waves in the complex lower atmosphere indicate that the 
energy in flu ctuations is distributed i ntermittently between spicules, loops, and the open- 
field corona ('De Pontieu et al."2007^; iTomczvk et"aDl2007l : iTomczvk and Mclntoshl l2009l : 
iFuiimura and Tsuneta 2009) . Magnetic flux tubes that thread the upper chromosphere and 
low corona support a wide r ange of possible non l inear couplings be tween compressible and 
incompressible modes (e.g.. lBogdan et alj|2003l : iHasan et al.ll2005h . For example, acoustic 
waves from non-magnetic regions of the chromosphere may encounter a swaying flux tube 
and be converted into Alfven or fast-mode waves as they pass through the tube. The study 
of energy transport across small-scale interfaces in the solar atmosphere may be crucial to 
producing more accurate and predictive models. 

Finally, it is important for future models to take a ccount of the k inetic and multi-fluid na- 
ture of coronal heating and solar wind acceleration JMarschlBOOq) . In the WTD paradigm, 
for example, a description of the large-scale energy flux injected into a turbulent cascade 
is only the first chapter in the story. Better descriptions of the anisotropic cascade pro- 
cess, the eventual kinetic dissipation of the fluctuations, and the subsequent energization 
of electrons, protons, and heavy ions are needed to complete the picture. Remote-sensing 
measurements of strong preferential heating and acceleration for heavy ions (e.g., 0+^) 
in coronal holes have spurred a great deal of the oretical work in this direction (see, e.g., 
iHoUweg and Isenberdl2002l : ICranmenl2002ll2009l) . A proper accounting of these kinetic ef- 
fects will lead to more concrete predictions for measurements to be made by space missions 
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such as Solar Probe JMcComas et alJl2007f) and Solar Orbiter JMarsden and FleckiBOOTh . 
as well as next-generation ultraviolet coronag raph spectro scopy that could follow up on the 
successes of the UVCS instrument on SOHO JKohl et alfc ooe. 2008,). 
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